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The Einstein's field equations of Friedmann-Robertson- Walker universes 
filled with a dissipative fluid described by both the truncated and non- 



^ ■ truncated causal transport equations are analyzed using techniques from dy- 

■ namical systems theory. The equations of state, as well as the phase space, are 

■ different from those used in the recent literature. In the de Sitter expansion 
On ■ both the hydrodynamic approximation and the non-thermalizing condition 

can be fulfilled simultaneously. For A = these expansions turn out to be 
bX). stable provided a certain parameter of the fluid is lower than 1/2. The more 



general case A > is studied in detail as well. 



I. INTRODUCTION 



Recently, isotropic spatially homogeneous viscous cosmological models have been inves- 
tigated using the causal (truncated and nontruncated) Israel-Stewart theory of irreversible 
processes, to modelize the bulk viscous transport^"!. It is known that dissipative pro- 
cesses may play a crucial role in the evolution of relativistic fluids both in cosmology and 
in high-energy astrophysical phenomena. The most oftenly used theory to describe such 
irreversible processes has been long since the first-order non-causal Eckart's theory! which 
however suffers from serious pathologies and drawbacks, i.e., superluminal velocities and 
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instabilitieJ-i. In the late sixties MiilleiS proposed a second order theory in which the 
entropy flow depended on the dissipative variables besides the equilibrium ones. Israel and 
Stewartil~i, and Pavon et a/.0 developed a fully relativistic formulation on that basis, the 
so-called extended or transient thermodynamics (see ref. 11 for a recent and comprehensive 
review of the state of the art). 

Shortly after Israel's paper appeared, Belinskii et al^ applied it to a viscous cosmological 
fluid using the so-called truncated version, in which some divergence terms in the transport 
equations were neglected. Most of the papers dealing with viscous and/or heat conduct- 
ing cosmological models make use of such a truncated transport equation without stating 
clearly what the implications of such a simplification may be. Recently, some effort has been 
invested in analyzing to what extent the neglecting of the divergence terms can be justified 
from a physical point of vie wi'i. As far as we know, Hiscock and Salmonson0 were the first 
to raise this point in the cosmological context. These authors stressed the key importance 
of the usually neglected divergence terms when obtaining viscosity-driven inflationary solu- 
tions. However, it is now clear that to get realistic solutions to the Einstein's field equations, 
the role played by the equations of state relating the different thermodynamic quantities is 
crucial. Hence the claim inlli applies only to a Boltzmann gaslii. In fact, the difficulty in 
using the extended transport equations lies mainly in the occurrence of some additional 
unknown coefficients, whose explicit expressions must be obtained from techniques other 
than those coming from thermodynamics, either kinetic or fluctuation theorylll, more than 
in their intrinsic complexity. 

Few exact solutions have been found to the Einstein's field equations with a non-perfec 
fluid described by extended thermodynamics0'lli (ET for short). However, they were ob- 
tained under severe restrictions on the values for the free parameters in the transport equa- 
tions. Obviously, any further attempt to get a deeper insight on the possible behavior of 
the solutions must rely on an approximate analysis of the equations. This paper applies 
qualitative analysis techniques to the study of causal viscous Friedmann-Robertson- Walker 
(FRW) models with and without a positive cosmological constant. It is organized as follows. 



In section II we state the basic equations governing the models and discuss the equations 
of state to be used. In section III we apply the truncated version of ET whereas in section 
IV a corresponding analysis is carried out using the full version. In both null and 

a positive cosmological constant are considered in turn. Section V explores their dynamical 
consequences. Finally, section VI summarizes the main conclusions of the paper. 



II. BASIC EQUATIONS 

We restrict ourselves to a FRW space-time filled with an bulk viscous fluid and a positive 
cosmological constant A. The stress-energy tensor is 

Tab= {p + P + 'n)UaUb+ {p+U- A)gab, (l) 

where Ua is the four velocity, p the energy density, p the equilibrium pressure, 11 the bulk 
viscous pressure. Einstein's field equations for the spatially fiat case (the only one we adress 
in this paper) are 

r, K A 

3^ 3' 

3{H + H') = ~ip + 3P,jf)+A, (2) 

where H = R/R is the Hubble factor, R{t) the cosmic scale factor of the Robertson- Walker 
metric, Pg// = p + 11 and k, = SttG/c^. An upper dot denotes differentiation respect to time 
t. We assume the fluid obeys equations of state of the form 

C = «p'", p=(7-l)p, r=^, (3) 

P 

where a is a positive constant, and 7 the adiabatic index lying in the range 1 < 7 < 2 as 
the sound velocity Vs/c = 7 — 1 in the fluid must be lower than the speed of light. r(> 0) 
is the relaxation time for transient bulk viscous effects, i.e. the time the system takes in 
going back to equilibrium once the divergence of the four-velocity has been switched off. 
The causal evolution equation for bulk viscous pressure can be cast into the form! 
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n + rri = -3CH-Ku (^3H + ^-^-^y (4) 

where 6 = for the truncated theory and 6 = 1 for the full one. Since a dissipative expansion 
is non-thermalizing, the relaxation time must exceed the expansion rate H~^. This leads to 

< H (5) 

which is a condition that reduces the interval of values of 7 for which the model holds. As 
we shall see this restriction may be violated in the truncated theory as well as in the full 
theory when an ideal gas equation of state is assumed. This conflict can be circumvented 
by resorting to the expression for the speed of the viscous signals a = v'^/c^ ~ C/'^'P? which 
roughly implie^'lli 

C 

r = — , < a < 1, (6) 
ap 

and using (H) instead of (0c). 

Most of the stability analysis below will be carried out for the de Sitter solutions, H = 
const. As the universe undergoes a de Sitter expansion it could be argued that the hy- 
drodynamic description (absolutely needed in our study if the results are to hold) might 
break down. In order the fluid approach remains valid the mean collision time tcoi must 
be less (in fact much less) than the expansion rate, i.e., tcoi < H^^. From kinetic theory 
one has tcoi = ^/na where n is the particle number density and a the cross section for 
collisions. In general cr(T) is an increasing function of the temperature whereas, for a FRW 
universe, n oc 1/B?{= e~^^"^ for a de Sitter universe). From equations (@a), (^) and (^ 
the non-thermalizing condition (^) and the condition for the hydrodynamic approximation 
imply 

p'iHot .o^ m-l 



noCT a \K, 

where rio is a positive but otherwise arbitrary integration constant. Later it will be shown 
that the second inequality in (H) can be fulfilled when suitable values for the arbitrary 



parameters are chosen. Moreover, the first inequahty may hold for sufficiently early times 
(when the inflation era supposedly took place). As the temperature remains constant during 
this period the cross section a can be taken approximately constantS . 

Recently in performing the qualitative analysis of imperfect fluid cosmological models 
(see for instance ref. 1, 2, 20, 21) dimensionless equations of state were used in terms of the 
dimensionless variables x and y defined as 

X = 3p/0^ y = m/e\ (8) 

The equations of state (eq.(3.4a,b) in ref. 20) 

p/Q^=p,x^, C/0 = Cox", (9) 

with 0(= 3H) the expansion factor, coincide with (^,b) only for / = m = 1/2. Furthermore, 
for the spatially fiat FRW metric with A = 0, the case we are interested in, we have x = 1. 
Then the bulk viscous coefficient ( varies as irrespective of m, which restricts (|^) to just 
one case: m = 1/2 in (^). In this paper we shall consider only the spatially flat case {k = 0) 
which allows us to take {H, H) as suitable dynamical variables in the phase space. In this 
case it appears to be more natural, especially when the above comments are taking into 
account, to adopt the oftenly used equations of state (||) rather than (|]) in order to be able 
to compare our results with those in the literature for m ^ 1/2. Consequently, all the fixed 
points to be analyzed will correspond to either de Sitter or static spacetimes {X = H = 0) 
being the former physically relevant in inflationary models. If one is interested in studying 
non-flat FRW models the variables {H, H) become no longer suitable. In such a case those 
proposed by Coley0 are more convenient. 

III. QUALITATIVE ANALYSIS USING THE TRUNCATED THEORY 

From equations (|^), (0) and the expression (D for r we find for the Hubble factor the 
equation 
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{?,H^ - A)a 



3H 



0. 



where S = a(3/K)™' ^. Equation (|TOD can be recast into the form 

H = P{H,X) 
X = QiH,X) 



(10) 



(11) 



where 



P{H,X)=X, 



Q{H,X) = -?>^HX~-[H g 



1— m 



(3i72 _ A)a 



X 



3i/ 



(12) 



(13) 



The quahtative analysis begins by hnearizing the system (|ll|) for small perturbations - where 
the linear theory holds. Then the Jacobian matrix 



( P p \ 



V 



with Ph = dP/dH, etc., 



(14) 



Qh Qx 

can be constructed. The elements of this matrix must be evaluated at the equilibrium points 
{hi,Xi) (de Sitter and static solutions) which are found by solving the system P{hi,Xi) = 
Q{hi,Xi) = 0. After diagonalizing L and obtaining its eigenvalues we can decide about the 
type of fixed points and their stability. 

The analysis of the system (|12]), (|l^) for the two cases with A = and A > will be 
carried out in turn. 



i. A = 

• m = 1/2 

We have the trivial fixed point 



(0,0), 
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which corresponds to an unstable static model. However, this case does not make sense as, 
by Einstein's equation (0a), p = 0. If 7 and 6 fulfill the restriction 

] = 3, (15) 

there exists an infinity of fixed points (/io,0), where ho denotes an arbitrary positive real 
constant. In this case the fixed points are parallel stable straight lines. 
• m^l/2 

In the intervals < m < 1/2, 1/2 < m < 2 there are two fixed points 

(0,0), (/io,0), (16) 

with 



1 — 2m 



ft„=(fj , (17) 

whereas for m > 2 there is only one fixed point {ho,0). The discussion for the point (0,0) 
mimics that for m = 1/2. 

Let us define the auxiliar parameter 

For m < | — Si the equilibrium point [ho, 0) is an asymptotically stable focus, for m = i — Si 
it is an asymptotically stable degenerate node, whereas for m > | — Si two cases arise. If 
I — Si < m < |, then the equilibrium point is an asymptotically stable node, whereas if 
m > 1/2 it is a unstable saddle point. 

In the paper by Pavon et sligthly different techniques were used to analyze the case 
m 7^ 1/2 and a = 1. Their relevant parameter was our Si with a = 1. Our results agree with 
those of the mentioned reference (see §3.1 of ref. 22) providing a more accurate classification 
of the stability points. 



ii. A > 
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As we shall see, there are two fixed points: (h^jO) and (/if,0). From ([T0| ) it follows the 
equation for the fixed points 

^ \ 1— m 



{3hi - A) 



h 



3hi 



0, 



(18) 



which must be solved for different values of m. However, (|T^) has an obvious solution 
independent of m, = a/A/3, which can be shown to correspond to a saddle point. This 



solution will be ruled out however since it would imply that the energy density vanishes 
identically. The other solution will be analyzed for m = 0, ^, 1, in turn. 
• m = 

Setting to zero the big square parenthesis in (^) and solving the resulting equation one 
obtains 
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We define 



A 



2Ei - r 72 (s^)2 ■ 

For A < Ao the fixed point {h^, 0) is an asymptotically stable node. If A = Aq the fixed 
point is an asymptotically stable degenerate node whereas for A > Aq it is an asymptotically 
stable focus. 

• m = 1 /2 
Now is given by 



A/3 



For < 952/^2 ^ (^1^,0) is an asymptotically stable node whereas for 9^2/72 > 1 the 
fixed point is an asymptotically stable focus for any A > 0. 
• m = 1 

Now hi = 7/35. For A < 72/3^2 -^^e have a saddle fixed point, whereas for < A < Ai, 

where 
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the fixed point is an asymptotically stable node. For A = Ai, it is an asymptotically stable 
degenerate node whereas for A > Ai, it is an asymptotically stable focus. 

IV. QUALITATIVE ANALYSIS USING THE FULL THEORY 

Actually, a proper study of viscous phenomena in the frame of ET requires the use of the 
full equation (Q) (i.e. 6=1). The physical implications of neglecting the second term of @ 
have been analyzed in detail in ref. 3, 13. The use of requires an explicit expression for 
the temperature T in terms of other variables such as p and/or n. So far (with the exception 
oE) the expression adopted for T has been a power-law 

T = (19) 

where r > and /? > are constants, which is the simplest way to guarantee a positive heat 
capacityii. However, we shall see that standard thermodynamic relations restrict the range 
of r. Calvao et a/.@ found a general equation for the evolution of temperature when two 
equations of state 

P = p{T, n), 

P = p{T,n), (20) 

are given. However, their equation was obtained in the context of matter creation where H 
is reinterpreted as a non-equilibrium pressure associated to particle production. The same 
equation has been carefully analyzed i n§, it reads 



T 



{dp/dT)n _^ H 



(21) 



idp/dT)^ T{dp/dT\^_ 
Obviously, when the equations of state (^) are known the evolution of T is no longer free 
but fixed by (^Tj). However, only in very few cases these equations are explicitly known, as 
for instance in the case of a radiation gas or an ideal gasil . Equation (|19D generalizes in a 
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simple way the Stefan-Boltzmann (r = 1/4) equation which holds for a radiation-dominated 
fluid in equilibrium. Thus, in that case we get that both equations of state p and p (when a 
7-law is used) have T as the only independent variable, i.e., dp{p)/dT = dp{p)/dT. When 
the integrability conditions for the Gibbs equation are considered an useful and interesting 
relation foUowH 

f dp\ p + P T ( dp\ 



dn J rp n n \ dT ^ 



(22) 



which, by virtue of (^b) and (|19D , yields 

7-1 



(^ < r < 1/2), (23) 

7 

i.e., r is no longer an independent parameter (we are indebted to Roy Maartens for pointing 
us out this restriction). It has been argued! that the inequality r < 1 is reasonable from a 
physical point of view, since ultrarelativistic and cold non-relativistic matter have r = 1/4 
and r ~ 2/3, respectively. 

However, an alternative equation can be used for T instead of (|19|). It is well-known that 



a relativistic ideal monoatomic gas is described by the two equations of state p = nT and 
p = 3nT + m?M, where m is the mass of the particles and M the zero-order moment of the 
Maxwell-Boltzmann distribution function (we use units = 1, fcs being the Boltzmann 
constant). We see that the 7-law (7 constant) is not compatible with the equations of state 
of a monoatomic gas in equilibrium except for radiation (m = 0). In that case we have 
n ocT^ and the two equations of state for p and p reduce to the Stefan-Boltzmann equation 
and the 7-law with 7 = 4/3. 

In the remainder of this section the full viscous transport equation will be analyzed 
resorting to the two expressions for the temperature mentioned above: a power-law given 
by (^9|) and an ideal gas equation for p together with the 7-law defining p, i.e. 

nT 

p = nT, p= -— ^. (24) 

Note that now both T and n are independent variables and only in the equilibrium limit the 
particle number density depends exclusively on the temperature, n = n{T) (see comments 
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above). It must be stressed that the Stefan-Boltzmann equation together with an ideal gas 
equation of state (with n (xT^) imphes 11 = 0. So we conclude that out of equilibrium we 
are forced to adopt one of the two possibilities: (i) a power-law for T with no dependence 
on n at all; (ii) an ideal gas equation of state for the pressure together a 7-law, with n an 
independent variable on the same footing as T. Both approaches will be considered in turn. 



A. Potential law for the temperature 

Using equations (0), (|^), ([T9D and (p3D the equation governing the evolution of the 
Hubble factor reduces to 



H 



i. A = 

• m = 1 /2 

As in Section III only the case with 7 and 5 fulfilling the restriction 

7 3 



(25) 



5 2 



(2-7), 



(26) 



is physically meaningful. In such instance there exists an infinity of stable fixed points 
(/ii,0), with hi an arbitrary positive real number. The phase portrait are parallel stable 
straight lines. 
• m ^ 1 /2 



There are two fixed points, (0, 0) and (/ii, 0), in the interval m G [0, |) 1J(|, 2) where 



hi 



3^(2 - 7) \~ 
, 27 J 



(27) 



For m > 2 only the {hi, 0) fixed point occurs, which we analyze next as nothing new arises 
about the point (0,0). 
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Let us define the parameter 

_ [^{2-a) + 2af ^ 
- 8a72(2-7) 

For m < I — S2 the equihbrium point is an attractor in the phase space (asymptotically 
stable focus). For ^ — S2 < m < ^ we have asymptotically stable nodes instead. Finally if 
m > 1/2, the fixed point is a saddle. 

ii. A > 



From eq. (|25|) it follows the equation for the fixed points 

1— m 



{3h} - A) 



A 



0, 



(28) 



where, as in the truncated case, only the solutions vanishing the big square parenthesis make 
sense from a physical point of view. Eq. (|28|) will be solved only for three different values of 
m. In this case we must take into account the constraint ([23|) . 

• m = 
The fixed point is (/i2 , 0) with 



4A 



Defining the two new parameters 



2(1 -i; 



i[l-f + fP-2(l-i) 



and 



A. 



27^%, 7,2 l + S^ 



we see that for A < A2 the fixed point is an asymptotically stable node, whereas for A > A2 
it is an asymptotically stable focus. 

• m = 1/2 
The fixed point is 
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A/3 



so 



3(2 - 7) 
Let us introduce 

For S > Sq, {h2, 0) is found to be an asymptotically stable node, however ii 6 = do it is an 
asymptotically stable degenerate node, whereas for 5 < So the fixed point is an asymptoti- 
cally stable focus. In the radiation case - i.e. 7 = 4/3- 5o < and the fixed point is a stable 
node. 

• m — 1 
Here 



m - i) 



Let us define the parameter 



_ 7^(1 + 2S,) 



For A > A3 the fixed point is an asymptotically stable focus, and for A = A3 an asymp- 
totically stable degenerate node. 

Finally when A < A3 we can distinguish two subcases. Defining 



A* ^ 



we have a saddle point for 

A< a;, 

and an asymptotically stable node for 

A* < A < A3. 
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B. Ideal gas equation for the temperature 



In this subsection we shall study the specific behavior of the equilibrium points mak- 
ing use of the state equations (p4D. As the number of particles is conserved n obeys the 
conservation equation 



h + 3Hn = 0, 



(29) 



which leads to n oc R . The expression for the temperature 



T = - 



K no 



3;' 



where no > is a constant, follows easily. Using @, (|a,b), 
evolution of the Hubble factor, the equation 



S\ 3j 2 



l—m 



(30) 

and (^) we get, for the 



-6 



3if2 _ A 



0. 



(31) 



i. A = 



We have the same fixed points as in the truncated theory (see equation (p^TD). 
For m = 1/2 the discussion parallels that of the truncated theory. After linearizing the 
system and introducing the parameter 

the following will be discussed. For m < | — S3 the eigenvalues are complex and the equi- 
librium point is an attractor (asymptotically stable focus). For m = | — S3 there is a 
bifurcation point which is an asymptotically stable degenerate node. For m > 1/2 one has 
a saddle point. Finally, if| — S3<m<l/2 the fixed point results an asymptotically stable 
node. 
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ii. A > 



Now the fixed points (/i^, 0) are again the same as in the full theory using a power law 
for the temperature. 

• m = 

For any A > the fixed point is an asymptotically stable focus. 

• m = 1 /2 
Defining 

^' = lf^^ 

we note that if 6 lies in the interval < 6 < 6i the fixed point is an asymptotically stable 
focus. If 6 = 6i it is an asymptotically stable degenerate node, and if 5i < 5 < 7/8 an 
asymptotically stable node for any A > 0. 

• m = 1 
Let us define 

A4 = ^(1 + 2S3). 

If A < 7/35^ then the fixed point is a saddle point, but if 7/35^ < A < A4 it is an asymp- 
totically stable node. For A = A4 it is an asymptotically stable degenerate node, and for 
A > A4 an asymptotically stable focus. 

Non-thermalizing condition for dissipative de Sitter expansion 

i. A = 

From (1^) and (0) one finds 

H^-^"" < 1 (32) 

For the truncated and full theory using an ideal gas equation for T this condition reduces 
to 7 > 3a by virtue of (0). On the other hand, as the velocity of the viscous pulses, as 
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well as the speed of sound, cannot exceeds the speed of light (1 < 7 < 2) we obtain the 
restrictions on 7 and a. If a lies in the range < a < |, the two mentioned conditions 
amount to 1 < 7 < 2; whereas if | < a < | these restrictions imply 3a < 7 < 2. Finally, if 
I < a < 1, no 7 can fulfill both conditions. 

For the full theory with a power law for temperature one obtains the restriction 7 > 7c 
where 

6a 

~ 3a + 2' 

7 must meet simultaneously two conditions: 1 < 7 < 2 and 7 > 7c. For < a < |, these 
restrictions imply 1 < 7 < 2 (since 7c < 1); whereas for | < a < 1 one has 1 < 7c < 7 < 2. 
So for a < 1 the full theory with a power law for temperature always holds. 

ii. A > 

Instead of (^2]) we now have 

when a positive cosmological constant is present. 

For the fixed point {h^, 0) (that of the truncated theory) the restrictions for 7 are the 
same that in the truncated case with A = 0, whereas for the full theory the restriction (p3D, 
when applied to the point (/i^,0), coincides with that of the full theory using a power law 
for the temperature with vanishing A. 

V. DYNAMICAL CONSEQUENCES 

In this section we study the dynamical implications of linearizing the equation for H. 
This linearitzation allows one to obtain an analytical solution for R{t) near the equilibrium 
points. The matrix L is given by (|Iip and the system of differential equations to solve is 
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(h 


\ 


/ 


Ph 


Px 


\ 


( 


k\ 






V 


Qh 


Qx 


J 


V 





(34) 



where h = H — hi being {hi,Xi) the fixed points. Equation (|3^ can be written as 

h~Qxh-QHh = 0, (35) 

where in our model Ph = and Px = 1- The corresponding characteristic equation reads 

Qx±^Qjc + ^QH 

Alt — 

2 

which coincides with the equation for the eigenvalues of L. We perturbe the system around 
the de Sitter solution for t = 0, i.e., H(t = 0) = hi + e(0) and take H{t = 0) = e(0), as 
initial condition. 



A. Saddle points and nodes 

In the neighborhood of these points the discriminant A is positive and the eigenvalues 
A± are real and different. For det{L) < one has a saddle point, and for det{L) > a node. 
The solution of (135|) is 



h = cie^+* + cse 



with 



Cl 



C2 



A_ - A^ 



e(0)-e(0)AH 



A_-A+ 

Upon integration one has for the scale factor 



R{t) oc e^'^ exp 



\J ^ A_ 



(36) 



which shows superinflationary expansion if initial conditions are taken such that ci, C2 
are positive. This type of evolution for R{t) on time has been obtained previously in a 
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different contextEZl. In this case the ffuid when submited to a small perturbation, departs 
from equilibrium, expands much more rapidly than the de Sitter's. In the case of nodes, 
and when A+ + A_ = Qx is positive (negative), the node will be unstable (stable). 



B. Attractors and repellors 



We study here the behavior of the scale factor near a sink (an asymptotically stable 
attractor) and a source (an asympotically unstable repellor). The solutions of the charac- 
teristic equation are complex quantities 



A- 



Qx±i\ \ A 



where A = Q\ + 'iQn- Then the solution of 



IS 



(t + Cs). 



h = Cic 2 sin 

The integration constants can be determined through the initial conditions. They read 

2 



C2 



tan 



A 



6(0) J| A 



2e(0) - e(0)g. 



and 



Cl 



e(0) 



sm -^^5—02 



Integrating the equation for h one follows 



R{t) oc e^'^ exp 



2ci 



Qx sin 



A 



it + C2) 



— \/ 1 A I cos 



A 



-(t + C2) 



(37) 



For Qx < we have an attractor (the scale factor undergoes an oscillatory approach to the 
de Sitter solution) and for Qx > it is a source, i.e., the scale factor deviates from the de 
Sitter solution. 
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C. Degenerate nodes 



In this case A = and A+ = A_ = A = The solution of (|35|) is 

h = cie^* + C2te^*. 
Because of the initial conditions the integration constants are 

ci = e(0), and Ca = e(0) - Ae(0). 
Integration of the expression for h leads to 



hence R{t) approaches to or separates from the de Sitter solution depending on the sign of 
A. The rate of evolution is faster than in the de Sitter case. 



The weak energy condition (WEC) states that TabWW^ > 0, where Tab is the energy- 
momentum tensor given by (|l]) and W"" a generic timelike vector. In our model this condition 



The dominant energy condition (DEC) imposes TaiW^W^ > and —T"'^Wa to be a non- 
spacehke vector which is equivalent to Too > iTa^lii. This conditions is fulfilled in our case 
only if 




(38) 



D. Energy conditions 



reduces to 



p + A > 0. 



(39) 



P<P + n<p + 2A. 



(40) 



Finally, the strong energy condition (SEC) requires that TabWW^ + lT^ > which amounts 



to 



p - 2A + 3p + sn > 0. 



(41) 



19 



These conditions can be rewritten in terms of H and H as 



WEC: 



SEC: 



DEC: 



H <0 and 3H^ + H>0 



As occurs in the standard inflationary scenarios the de Sitter solutions with A > satisfy 



We have carried out a detailed analysis on the stability of de Sitter and static cosmolog- 
ical models, both in the truncated and full theory for the viscous transport equation (with 
and without a cosmological constant) . We have shown that the conditions for the hydrody- 
namic approach and the non-thermalizing condition in a de Sitter expansion can be fulfilled 
simultaneously for sufficiently early times. When no cosmological constant is considered 
the stability analysis for the de Sitter solutions leads to similar results in all the cases, i.e., 
the models are unstable only for m > 1/2. It is remarkable that this result holds for both 
the truncated and the full version of ET. On the other hand, when a positive cosmological 
constant is included we see that the models can be stable for m = 1 if A is bounded from 
below. It remains to be proved that this result holds for a generic m > 1/2 other than 1. 

We have stressed the fact that for a radiation gas a different thermodynamic approach 
exists depending on whether the particle number density is taken as an independent variable 
or not. In the case of a power-law for the temperature a relationship between 7 and r exists. 
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